******************
**** Figure 5 ****
******************

* Run table b7 to get estimates E *
do "$dofiles/table b7" 

* plot
coefplot E, keep(treatmentalldummy2 treatmentalldummy3 treatmentalldummy4 treatmentalldummy5 treatmentalldummy6 treatmentalldummy7) vertical yline(0, lcolor(black)) recast(bar) barwidth(0.9) lcolor(gs0) fcolor(*.5) color(gs6) ///
ciopts(recast(rcap)) citop  coeflabel(treatmentalldummy2="P-Bas" treatmentalldummy3="N-Bas" treatmentalldummy4="P-Req" treatmentalldummy5="P-Attr" treatmentalldummy6="N-Req" treatmentalldummy7="N-Attr") format(%9.2f) ///
msymbol(t) levels(95) graphregion(color(white)) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) ///
ylabel(-0.15  -0.1 -0.05  0  0.05  0.1  0.15  0.2 " ", angle(0)) ///
ytitle(ITT Effect on Collection Efficiency) title(Post-Intervention Treatment Effects, color(black)) subtitle("September 2017 - May 2018") ///
groups(treatmentalldummy2 & treatmentalldummy3  = "{bf:BASIC}" treatmentalldummy4 & treatmentalldummy5 = "{bf:PIDA}" treatmentalldummy6 & treatmentalldummy7 = "{bf:NIDA}")

graph export "$filepath/output_logs/figure 5.pdf", replace

